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This paper presents an exploration of two fundamental scaling relations of 
spiral galaxies, the luminosity-rotation speed (or Tully-Fisher; TF) relation, and 
the size-luminosity (SL) relation, and the dependences of their scatter, at red and 
infrared bands. We verify that the observed virial relations of disk galaxies are 
given by Vobs oc Lf 31 and Rd oc Lf 33 using distance-redshift surveys of high sur- 
face brightness (HSB) and low surface brightness (LSB) non-interacting galaxies. 
These results imply that the galaxy surface brightness £ oc R d oc V ohs oc M*/ 3 . 
The collected surveys provide accurate /-band luminosities, disk scale lengths, 
circular velocities, and, in some cases, color. Various issues regarding the scatter 
of scaling relations in blue to near-infrared bands are being re-examined with 
accurate JHK Kron luminosities, effective radii, and colors from the 2MASS 
database. We derive the first extensive J-band TF and SL relations. At a given 
infrared luminosity, the TF velocity residuals are correlated with infrared color, 
which in turn is determined by the variations in galaxy formation ages and dark 
halo concentrations; these residuals are fully independent of surface brightness 
and other tested galaxy observables. We verify that TF relations (TFRs) for 
various morphological types have different zero points, but a common slope, such 
that early-type (redder) disks rotate faster than later-type (bluer) systems of 
the same luminosity. The morphological type dependence of the TFR is a di- 
rect consequence of the more fundamental scatter dependence on color, which 
itself is related to the star formation history of a galaxy. The scatter of the 
SL relation is mostly dominated by surface brightness, but color also plays a 
small role. The observed systematic variations of disk size with color, for a 
given luminosity, are weaker than expected, perhaps as a result of uncertainty 
in disk scale length measurements. The TF and SL residuals for HSB and LSB 
galaxies are weakly correlated with d\ogV(L) / d\og R(L) = —0.07 ± 0.05, in 
agreement with the earlier claim by Courteau & Rix (1999, ApJ, 513, 561), and 
d\ogR(V) I dXogLiy) = 0.12 ± 0.03. The former result suggests that spiral 
disks of all color, brightness and barredness may be, on average, dominated by 
dark matter even in their inner parts (the so-called "sub-maximal disk" solution 
at R > 2.2 disk scale lengths). The observed scaling relations of disk galaxies are 
likely the result of the simplest scenario of galaxy formation in the standard cos- 
mological picture, but their color dependence remains a challenge to hierarchical 
structure formation models. These relations and their derivatives are analysed in 
terms of, and serve as stringent constraints for, galaxy formation models in our 
companion paper. 

Subject headings: galaxies: dynamics — galaxies: formation — galaxies: kinemat- 
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ics — galaxies: spirals — galaxies: structure — dark matter 

1. Introduction 

Understanding the origin and nature of galaxy scaling relations is a fundamental quest 
of any successful theory of galaxy formation. The success of a particular theory will be 
judged by its ability to predict the slope, scatter, and zero-point of any robust galaxy scaling 
relation at any particular wavelength. Some observed scaling relations in spiral galaxies, 
based on their size, luminosity, and rotation speed, can be reproduced individually to fairly 
good accuracy by invoking galaxy formation models that include virial equilibrium after 
dissipational collapse of spherical cold dark matter (CDM) halos and angular momentum 
conservation (e.g. Mo, Mao, & White 1998, hereafter MMW98; van den Bosch 1998, 2000, 
hereafter collectively as vdBOO; Navarro & Steinmetz 2000, hereafter NS00; Firmani & Avila- 
Reese 2000; hereafter FAR00). 

One of the most firmly established empirical scaling relation of disk galaxies is the Tully- 
Fisher relation (TFR; Tully & Fisher 1977); a tight correlation between the total luminosity 
and the rotation speed of a disk galaxy. However, to date, no single CDM-based model 
of galaxy formation can simultaneously reproduce the slope, zero-point, scatter, and color 
trends of the TFR, match the shape and normalization of the luminosity function, and 
explain the sizes, colors, and metallicity of disk galaxies (see, e.g,. vdBOO; Bell et al. 2003). 
In addition, simultaneously accounting for the mass and angular momentum distribution 
of spiral galaxies in a gas dynamical context remains a major challenge for hierarchical 
formation models (Navarro & White 1994; Bullock et al. 2001b; van den Bosch et al. 2002b). 
A complete theory of galaxy scaling relations awaits a fuller understanding of structure 
forming mechanisms and evolutionary processes (e.g. star formation, merging, feedback, 
and cooling prescriptions) in galaxies (see Somerville & Primack 1999 for a comparison of 
"recipes" among competing theories.) Likewise, the fine-tuning of these galaxy formation 
and evolutionary models demands a careful examination of empirical scaling relations of 
galaxies. 

In order to set up a framework for the study of galaxy scaling relations, we examine 
the correlations of parameters related by the virial theorem, V 2 oc Mj R, for bright galaxies 
(i.e. where feedback effects are minimal). We consider three fundamental observables for 
each disk galaxy: the total luminosity L, the stellar scale length i?d of the exponential disk, 
and the observed circular velocity Vobs- The stellar mass, Md, can be estimated from the 
luminosity by assuming a stellar mass-to-light ratio, Yd = M^/L. The size-luminosity (SL) 
relation of galaxy disks is also expressed as L oc where S is the surface brightness. 
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Key to mapping fundamental dynamical trends in spiral galaxies, the measurement of 
TF and SL relations and detection of their correlated residuals require velocity amplitudes 
measured at a suitably chosen radius representative of the flat part of resolved rotation 
curves, red/infrared magnitudes to minimize extinction and population effects, accurate 
disk scale lengths, and, ideally, color terms with a broad baseline (e.g. B — K) to test for 
Yd variations in the stellar population and extinction effects. Our goal in this paper is thus 
to assemble such a data base. 

A study of scaling relations in irregular and spiral galaxies by Salpeter & Hoffman (1996; 
hereafter SH96) yielded the correlations Lb oc R 2 - 68 oc \/ 3 bs 73 oc M^ 35 oc M]^ 6 , where M H , 
Mdy n , and R, are the Hi and dynamical masses, and a characteristic radius, respectively. 
The blue luminosities, as used in that study, are notoriously sensitive to dust extinction 
and stellar population effects and dynamical effects cannot be simply isolated. A new study 
of scaling relations in the (near-)infrared would provide more robust dynamical constraints 
to galaxy formation models. While the TFR has been examined at nearly all optical-IR 
wavelengths (e.g. Strauss & Willick 1995; Verheijen 2001, hereafter V01), comparatively few 
multi-wavelength analyses of the SL relation (SLR) of spiral galaxies have been done so far 
(SH96; MMW98; Shen et al. 2003). This is partly because the accurate disk scale lengths 
needed to calibrate the SLR, at any wavelength, have only recently become available for large 
databases (e.g. Courteau 1996; Dale et al. 1999; MacArthur et al. 2003, hereafter MCH03). 
The recent availability of light profile decompositions for very large galaxy databases (e.g. 
Sloan Digital Sky Survey and Two Micron All-Sky Survey; §2) heralds a new era for the study 
of galaxy scaling relations (e.g. Shen et al. 2003). We will return to considerations about 
the SLR in §2.1. The remainder of this section focuses on the TFR and the dependences of 
its scatter on galaxy observables. 

The empirical TFR is expressed as 

L oc KTbs (1) 

with the near-IR log-slope a ~ 3.0 ±0.4 (Willick et al. 1997; Giovanelli et al. 1997, hereafter 
G97; Courteau et al. 2000; V01). Reported values of the log-slope a range from 2.8 in the blue 
to 4.0 in the infrared (Willick et al. 1997; Tully & Pierce 2000; V01) for both high and low 
surface brightness galaxies (Zwaan et al. 1995; V01). Log-slopes steeper than a ~ 3.5 in the 
infrared typically result from small samples and excessive pruning on the basis of idealized 
morphology or kinematics, a narrow range of inclinations, redshift cutoffs, etc. (Bernstein 
et al. 1994; V01; Kannappan, Fabricant, & Franx 2002, hereafter KFF02). The slope, scatter, 
and zero-point of blue TFRs are predominantly dominated by stellar population and dust 
extinction effects (e.g. Aaronson & Mould 1983; C97; G97; Willick et al. 1997; Tully & Pierce 
2000) and on the techniques used to recover the major observables and fitting for fundamental 
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relations (e.g. Strauss & Willick 1995; C97; V01; Bell & de Jong 2001; KFF02). Because 
we are mainly interested in masses, rather than luminosities, we do not concern ourselves 
with TFRs and other scaling relations measured at blue wavelengths. We show in Appendix 
A that the fundamental form of the TFR, based on dynamical principles only, is given by 
a = 3, that is L oc V^, s . 

The modern interpretation of the TFR is that of a correlation between the total baryonic 
mass of a galaxy, inferred via its infrared luminosity and a stellar mass-to-light ratio and 
total gas mass (Hi + He), and its total mass inferred from the asymptotic circular velocity 
of the galaxy disk (McGaugh et al. 2000; Bell & de Jong 2001; V01). The "baryonic" TFR 
is expressed as 

M hary oc V2r (2) 

Since disk gas mass fractions typically increase with decreasing luminosity (i.e., Mc- 
Gaugh & de Blok 1997), one typically has that ab ary < a. 

While the log-slope of the TFR can be reproduced fairly well by most CDM-based 
structure formation models (e.g. MMW98; vdBOO; NS00), the predicted scatter can be large 
compared to the inferred "cosmic" scatter of 0.25 mag in red/infrared bands (Willick 
et al. 1996; V01) and interpretations about its dependence differ (see below). 

Besides the basic understanding of the slopes of galaxy scaling relations, the dependence 
of their scatter has also been addressed by many, especially for the TFR (e.g. Aaronson & 
Mould 1983; Giraud 1986; Rhee 1996; Willick et al. 1997; KFF02), and can be used to set 
stringent constraints on structure formation models (Courteau & Rix 1999, hereafter CR99; 
Heavens & Jimenez 1999; FAR00; NS00; V01; Buchalter, Jimenez, & Kamionkowski 2001; 
Shen, Mo, & Shu 2002). 

While various trends in the scatter of the blue TFR have been reported in the past, no 
correlations of the infrared TF residuals with inclination, size, concentration, gas fraction, 
or far infrared luminosity have thus far been reported (Aaronson & Mould 1983; V01). The 
dependence of near-IR TFR scatter on color and surface brightness is however still a matter 
of contention that we discuss in §2.2. The SLR scatter is also addressed in §2.2. 

The study of scaling relations in galaxies has benefited from the two-pronged applica- 
tion of the TFR for the purposes of: (i) Estimating relative distances to measure deviations 
from the mean Hubble flow (see, e.g., Strauss & Willick 1995 and the reviews in the "Cosmic 
Flows 1999" proceedings by Courteau, Strauss, & Willick 1999); and (ii) testing galaxy for- 
mation and evolution models (Dalcanton, Spergel & Summers 1997; MMW98; vdBOO). The 
philosophy of sample selection and calibration differs in both cases. For cosmic flow analyses, 
the calibration and science samples must be pruned mostly on the basis of morphology and 
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visual appearance in order to minimize systematic errors and thus ensure the smallest possi- 
ble magnitude (distance) error. Since TF scatter depends strongly on the slope of the TF, it 
is found that the combination of steepness, magnitude errors, extinction correction, and sky 
stability favors red (RSzI) bands for smallest distance errors and cosmic flows applications 
(Courteau 1997, hereafter C97; V01). The accuracy of bulk flow solutions also depends on 
the size of the sample. In order to collect large enough samples, TF calibrations for flow 
studies rely mostly on Hi line widths or Ha rotation curves that can be collected relatively 
quickly on modest aperture telescopes. These rotation measures typically sample the disk 
rotation out to 2 to 3 disk scale lengths (C97). 

By contrast, use of the TFR as a test bed for galaxy formation models requires the widest 
range of morphological types, to sample all structural properties, and that extinction and 
stellar population effects be minimized to isolate genuine dynamical correlations. The nearly 
dust-insensitive i^-band is thus the one of choice for such applications (V01) 1 . Rotation 
velocities are preferably extracted from fully resolved Hi rotation curves obtained using 
aperture synthesis maps that sample the disk rotation out to 4 to 5 disk scale lengths. 

Ideally, the study of scaling relations should rely on homogeneous samples assembled 
with the very purpose of testing for broad structural and dynamical differences; at the mo- 
ment, we must contend with the more finely pruned heterogeneous samples of late-type 
spirals that have been collected during the last decade mostly for flow studies. These data, 
which include near-infrared luminosities, and sometimes colors, for large samples of galax- 
ies, still enable us to characterize the dependence of scaling relations and examine various 
constraints of structure formation models, in ways hitherto unsettled. 

In Courteau et al. (2003), we showed that barred and unbarred galaxies have similar 
physical properties and that they share the same TFR. As a natural extension of this study 
and CR99, in §2.2 we use the extensive all-sky distance-redshift catalogs presented in §2 
to characterize the scatter of the TF and SL relations in terms of galaxy observables for a 
broad suite of galaxy types. In §3 we bolster the notions of surface brightness independence 
and color dependence of the TFR, as well as that of surface brightness dependence for the 
SLR, for disk galaxies. This prompts a renewed examination of the correlated scatter study 
of CR99 which favored high dark-to-luminous mass fractions in galaxy disk's interiors. 

In §4 we develop a simple analytical interpretation for the origin of the galaxy virial 
relations, such as the TF and SL relations, in a cosmological setting. The analysis of the 
scatter in the scaling relations, and a discussion about the implications of our results in terms 



*We find in §2.3 a slightly tighter 2MASS TFR at J-band rather than at K; the ff-band TFR is hardest 
to control due to airglow fluctuations (Jarrett et al. 2003). 
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of existing galaxy formation models, are presented in our companion paper (S. Courteau 
et al. 2004a, in preparation; hereafter Paper II). 

2. Available Samples 

We consider four samples for which accurate galaxy observables, including crucial rota- 
tional velocities, are available. These include (a) the large /-band survey of galaxy distances 
for bright field spirals in the Southern sky by Mathewson, Ford, & Buchhorn (1992; here- 
after "MAT"); the all-sky /-band TF surveys by (b) Dale et al. (1999; hereafter "SCII") 
and (c) Courteau et al. (2000; 2004b, in preparation; hereafter "Shellflow"); and (d) the 
multi-band BVRK TF survey of Ursa Major cluster galaxies by Tully et al. (1996) and 
V01 (hereafter "UMa"). The first three TF surveys were originally designed to map the 
convergence of the velocity field on ~ 60/i -1 Mpc scales, with the SCII and Shellflow studies 
paying special attention to TF calibration errors between different observatories. By design, 
these surveys favor late- type galaxies with inclinations on the sky greater than 32°; most 
have i ~ 60°. 

Properties for each sample are given in Table 1. These include in Col. (2), the number of 
galaxies from the original database that have a full complement of useful observables; Col. (3), 
the nature of the sampled galaxies (cluster or field; the predominantly "field" surveys of MAT 
and Shellflow include a small fraction of cluster galaxies); Col. (4), the digital photometric 
coverage. -B-band magnitudes for all the MAT galaxies were extracted from the RC3 (de 
Vaucouleurs et al. 1991). We have Sloan Digital Sky Survey (York et al. 2000, hereafter 
SDSS) g — % colors for 39 SCII galaxies (see Appendix B) and JHK Kron magnitudes and 
colors from the Two Micron All-Sky Survey (Skrutskie et al. 1997; hereafter 2MASS) for 
the brightest ~ 400 SCII galaxies; Col. (5) the magnitude or diameter limits of the original 
catalog; Col. (6), the redshift limits of the survey in kms -1 . For UMa, we list the adopted 
cluster distance in Mpc; and in Col. (7), the rotation measure, either extracted from Hi line 
mapping or Ha spatially resolved rotation curves. While MAT gives both Hi and Ha rotation 
measures we use only the sub-sample of galaxies with available Hi , which is largest. Our 
results for the MAT sample do not depend upon this particular choice of rotation measures. 
For SCII, rotation velocities were obtained using either Ha long-slit spectra or Hi line profiles. 
These two rotation measures are comparable as long as the rotation curve is flat beyond 2 
disk scale lengths and the Hi surface density drops rapidly in the outer disk (C97). Rotation 
velocities for the Shellflow sample were measured at R ~ 2.2/?d from resolved Ha rotation 
curves. For the UMa database, we use the "Hi " sample of 44 galaxies with Hi line widths 
(V01). 
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Disk scale length measurements were computed in somewhat different ways for each 
sample. For the MAT galaxies, we have used the two-dimensional bulge-disk (B/D) de- 
compositions of MAT /-band images by Byun & Freeman (1995). These decompositions 
assumed a de Vaucouleurs bulge profile which is not ideal for late-type galaxies (MCH03) 
and therefore the disk central surface brightnesses are likely biased low. The disk scale 
lengths are less affected, in a relative sense due to the broader baseline of the disk modeled 
as an exponential profile, and deemed adequate for the purpose of our study (see §2.1). 
SCII disk scale lengths were obtained by fitting a straight line to the exponential part of 
the /-band surface brightness profile from ~21 /-mag arcsec~ 2 to ~25 /-mag arcsec -2 (the 
so-called "marking the disk" technique) and corrected for projection effects according to 
Rd = Z?d[l + 0.41og(a/6)], where a and b are the semi-major and semi-minor axes of the 
disk (Dale et al. 1999; G97). The correction for projection is small and somewhat uncertain 
but the choice of observed or deprojected scale lengths for the SCII galaxies does not affect 
our final results. Shellflow scale lengths were extracted from one-dimensional B/D decom- 
positions of azimuthally-averaged /-band surface brightness profiles (MCH03; S. Courteau 
et al. 2004; in preparation). These fits account for a Sersic bulge and an exponential disk. 
Disk scale lengths for the UMa sample were measured from a "marking the disk" technique 
but the fit baseline is unspecified and erratic fits are reported in Tully et al. (1996). The 
latter could be due to inclusion of bulge light in the disk fit (scale lengths biased low) for 
the brighter galaxies and sky domination (scale lengths biased high) for the fainter galaxies. 
The latter signature is detected when we compare the scale length measurements for the 
four samples in §2.1. 

Physical parameters are computed using Ho=72 km sec" 1 Mpc -1 . This matches the 
revised distance to the UMa cluster from 15.5 Mpc to 18.6 Mpc, based on new Cepheid 
distances from the HST Key Project (Tully & Pierce 2000). The exact choice of distance 
scale does not affect our conclusions so long as it is consistent for all samples in order to put 
absolute scales on equal footing (luminosities in L & and scale lengths in kpc). 

Hubble types for all but the SCII galaxies were obtained from the heterogeneous NASA 
Extragalactic Data (NED) Base. Most SCII galaxies were not classified in NED and their 
morphological types were determined in a homogeneous fashion via a combination of eye- 
ball examination, B/D ratio and/or concentration index (Dale et al. 1999). Interacting and 
disturbed galaxies were rejected in all samples. Together, the MAT, SCII, Shellflow, and 
UMa samples combine for a total of 1,750 separate entries; some repeat measurements within 
a given sample exist but all observations are independent. 
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2.1. Mean Parameter Relations 

We consider the projected distribution in each of the planes defined by a pair of the three 
log virial variables, log Rd, log V b s , and log L; these are shown in Figs. 1-4 for the SCII, MAT, 
Shellflow, and UMa galaxies. The axis limits are the same in all figures for trend matching 
between each sample. Different symbols identify the full range of spiral Hubble types, as a 
proxy for accurate optical colors which are not available for the full MAT and SCII samples. 
The solid lines correspond to one dimensional linear regressions to (a) log V Q b s on log Lj, (b) 
logi?d on log Lj, and (c) logi?d on logV^bs, given the measurement errors in the first variable 
in each case when available (or otherwise using unweighted fits). The fits for the SCII sample 
are used as fiducials against the other, less complete, samples and are shown as a dotted 
line on the virial projections for each sample. For each relation, we fit the functional form 
y(x) — y + a y (x — x), where y and x stand for any combination of log V b s , log Rd, and log L; 
the upper bar denotes the median value of each variable. To achieve a robust fit to the data 
[y(i),x(i)], we varied a y to minimize the data— model deviation (rather than the squared 
difference). We repeat the fits for all the log quantities keeping either the total luminosity, 
stellar scale length, or maximum velocity fixed. This fitting technique follows the treament 
in CR99. We have also computed standard least-squares fits which yield similar results. 
Log-slopes and the Pearson linear correlation coefficients r for each parameter combination 
are reported in Table 2 (we do not concern ourselves with calibration zero-point differences 
in this Paper). 

The regressions were performed using the full range of luminosities (when the regression 
was on logL) and velocities (when the regression was on logV^s)- It has been proposed 
that the scaling relations may be different for the fainter, LSB galaxies (e.g., Kauffmann 
et al. 2003; Shen et al. 2003), where supernova feedback effects are likely to be important 
(e.g., Dekel & Silk 1986; Dekel & Woo 2003). Such a departure is, however, not apparent in 
the samples (SCII and MAT) that most closely probe the LSB regime (see Figs. 1 & 2). 

The agreement in Table 2 between the log-slopes for different samples is relatively good, 
though differences exist. All samples are well matched by a TFR, V oc L a , with log-slope 
ai = 0.31. Our TF fits are also a close match to those published by the original authors (as 
M a bs oc logVobs), even though the minimization techniques can be quite different. The Pear- 
son correlation coefficient of the distribution (not considering the errors) averages r = 0.90 
for the TFR and the conditional distribution of V ohs at a given L is roughly log-normal. The 
correlations for the SL and size-velocity (SV) relations are weaker, averaging r ~ 0.6. The 
size distributions are less robust than the TFR owing to a number of factors such as, (a) small 
range of galaxy sizes in some samples (Shellflow, UMa), (b) lack of a uniform and universal 
definition of scale length (MCH03), and (c) intrinsic scatter due to the natural dispersion of 
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the spin parameter, A (§4). Whether scale lengths are measured from B/D decompositions 
or "marked" over a specified range of surface brightnesses (while omitting the bulge region 
or not) can yield scale differences greater than 20% (MCH03). The overlap between the four 
samples is too small (in some cases non-existant) to calibrate such systematic errors. Thus 
it is difficult to determine how much of the lower correlation coefficients in the SL and SV 
relations is due to genuine dynamical processes and/or inadequate scale lengths. In spite of 
potential pathologies with the MAT scale lengths (see §2), the scalings for this sample and 
SC11 are comfortably close. For the sake of uniformity, we will base our final analysis of 
the virial relations of disk galaxies on the SCII sample alone, being the most extensive and 
robust sample considered here. Furthermore, unlike our results for some of the other, less 
complete and thus noisier, data sets our fits for SCII are self-consistent, within the errors. 
If V oc L a and R oc L 13 and V oc i? 7 , then self-consistency requires that a = f3^. This is 
indeed the case for SCII, our most complete sample. 

In summary, we find the following mean power- law scaling relations at J-band for spiral 
galaxies: 

KbsOcL?- 31 , R d oc L°j U , i? d ocKbs 00 - (3) 

The uncertainties in the log-slopes for the SCII sample, based on measurement errors, are 
±0.02, ±0.02, and ±0.09, respectively. 

Fig. 5 shows the combined TF and SL relations for the four samples. This operation 
formally requires that the magnitudes, rotation velocities, and scale lengths be perfectly 
homogeneous. However, zero-pointing differences exist between the various samples thus 
preventing a direct merging of all samples into one large homogeneous catalog for our study 
of scaling relations. The rest of our investigation thus relies on the separate examination of 
each data set. 

2.2. The Morphological/Color Dependence of Galaxy Scaling Relations 

The morphological type dependence of the TFR is highlighted in Figs. 1-4 by different 
point types and colors. This dependence, first noted by Roberts (1978), and revisited by 
Aaronson & Mould (1983), Rubin et al. (1985), Giraud (1986), Pierce & Tully (1988) and 
G97 - to name a few - goes in the general sense of early-type galaxies rotating faster than 
later-types at a given optical luminosity. This trend is confirmed for each sample in the top 
panels of Figs. 1-4 where the early-type and late-type spirals, represented by red stellated 
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points and black filled triangles, are located above and below the fiducial fits, respectively 2 . 
Recent optical TF relations of nearby spirals (e.g. KFF02) and of distant spirals (B. Weiner 
et al. 2003 [DEEP project], in preparation) covering a wide range of morphologies support 
the Hubble-type dependence of the optical TFR, up to at least z ~ 1. 

Figs. 6 & 7 show histograms of the luminosity residuals from the mean TF relation as 
a function of Hubble type for the Shellflow and SCII samples (similar results are obtained 
for MAT and UMa). By selection, the mean TF fit is dominated by Sb-Sc galaxies. For 
the Shellflow sample, early and late-type spirals differ by 0.08 dex and -0.25 dex from the 
mean TFR respectively. This is accentuated in the SCII sample (Fig. 7) with TF luminosity 
residuals of +0.17 dex and -0.33 dex for early and late-types, respectively. Dale et al. (1999) 
applied a morphological-type correction to their galaxy magnitudes for the construction of 
their TF template. Their corrections to total /-band magnitudes for T<Sab, Sb, and T>Sbc 
are -0.27, -0.11, and 0.00 mag, respectively, which is somewhat larger than the ones we find 
for their data. 

While Hubble-type classification depends on the bandpass of selection, a morphological- 
type dependence of the TFR, even at near-infrared bands, is undeniable. The morphology 
dependence of the SV and SL relations is weaker, as one would expect based on the rela- 
tively larger errors of disk scale lengths compared to circular velocity measurements. The 
distributions of scale lengths in the SL diagrams show a slight excess for bluer (later type) 
galaxies over the mean at a given luminosity, suggesting that bluer galaxies are bigger than 
the mean. Prejudice from theoretical models (e.g. FAR00; vdBOO) and our own simple 
scenario of disk formation (§4) would have us favor a scenario where redder objects are both 
faster rotators and more compact, but observational evidence for the latter is weak. This 
could be the result of inhomogeneous disk scale length measurements, and especially for the 
fainter systems for which sky errors are more conspicuous. The slight excess of disk scale 
lengths over the mean at the faint end of the SLR could be understood if sky subraction 
was systematically under-estimated. In that respect alone, an extensive program to measure 
accurate disk scale lengths and half-light radii for a large collection of spiral galaxies based 
on one unique reduction method and high signal-to-noise data is badly needed. 

Because morphological classification is a subjective measure, we look at more objective 
observables to determine what physical parameters drive the TFR scatter. Figs. 8 and 9 
display the TFR scatter in magnitude against five galaxy observables for the four samples. 



2 At the lenticular end of the spiral sequence, SO galaxies are consistent with our qualitative findings with 
a marked offset from the mean TFR (of Pierce & Tully 1992) but these galaxies also show a steeper TFR 
than early-type spirals (Neistein et al. 1999; Hinz, Rix, & Bernstein 2001; Mathieu, Merrifield, & Kuijken 
2002). The steeper TFRs may be a consequence of sample selection and/or just the nature of SO galaxies. 
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Besides the RC3 and SDSS colors extracted for the MAT and SCII samples respectively 
(see Appendix B for SCII), the data were all collected from the original literature. The 
concentration indices (CIs) for Shellflow and SCII galaxies correspond to the ratio of radii 
containing 75 and 25 percent of the total light (CI72). CIs for UMa galaxies were computed 
at radii containing 80 and 20 percent of total light (CI82; Tully et al. 1996). Note that all 
exponential disks without a bulge have the same value of CI72=2.8 and CI82=3.6, regardless 
of their total mass or scale length. For the MAT sample, we substituted missing CI data 
with B/D ratios. 

Looking at Figs. 8 & 9, we see that the TFR scatter (in magnitudes) is clearly indepen- 
dent of central surface brightness, /io, physical disk scale length, R^, and concentration index. 
Color and morphological types, on the other hand, are both correlated with TF magnitude 
residuals. Given that color and morphological type are themselves strongly correlated (see 
Fig. 10), this is consistent with a picture in which color is the fundamental driver of scatter 
in the TF relation. Fig. 10 shows that the tightest correlation with morphological type is 
indeed with color. Fig. 10 also shows that velocity depends more strongly upon morpholog- 
ical type (color) than disk size, an important observation that we return to in §3 and Paper 
II. The correlation of morphological types with concentration index is poorest, highlighting 
the relative inadequacy of the CI test to discriminate between different spiral Hubble types. 

According to V01, all correlations of TF residuals with galaxy observables detected in 
blue bands, including color, vanish in the infrared. This important conclusion however rests 
on the study of a small data sample (less than 44 galaxies) and deserves closer attention. 
Contrary to findings in Aaronson & Mould (1983) and V01, Rubin et al. (1985) found 
that the morphological type dependence of the TFR was reduced, but not eliminated, at 
infrared H luminosities 3 . Thus color is potentially another fundamental parameter of galaxy 
structure and formation, in addition to size, luminosity and rotation speed. Until recently, 
the CCD luminosity measurements in large galaxy samples with available line widths were 



3 Unlike V01 and Aaronson & Mould (1983), Rubin et al. (1985) found a similar morphology (color) 
dependence of the TFR, even at _ff-band. Such a dependence could be an artifact of rotational speed 
measurements with optical spectra. The shape of optical rotation curves is often dominated by the inner 
rising part of an early-type spiral or the slowly rising part of a later-type spiral of the same luminosity. 
However, the presence or absence of a bulge or bar has a weaker effect on the overall extended shape of a 
resolved HI rotation curve. Thus the asymptotic rotational speed measured from HI synthesis maps or Ha 
velocity fields from say, 2 disk scale lengths to the last detected velocity point (with acceptable signal-to- 
noise), will not be subject to the vagaries of the inner rotation curve in the same way as optical rotation 
curve measurements (e.g. C97, V01). So long as velocity measurements are made beyond R = 2.2i?<j, the 
rotation speed should be unaffected by a bulge or bar (CR99; Courteau et al. 2003), and the trends detected 
by Rubin et al. and ourselves are robust. 
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only collected in a single band (typically r or I), thus thwarting any systematic test for 
color dependence. Recent scatter studies based on the small, pruned, UMa galaxy sample 
(Heavens & Jimenez 1999; V01) can now be duplicated for much larger samples, as we do 
here. 

The recent i?-band TF study of the Nearby Field Galaxy Survey (Jansen et al. 2000) 
by KFF02 identified B — R color and Ha equivalent width as the main drivers of the TFR 
scatter. Color and Ha equivalent width are both tributary of star formation histories though 
the former depends both on the integrated and instantaneous star formation while the latter 
is a function of the current star formation rate alone. Taken together with our results, it 
thus seems that star formation history is the fundamental driver of the TF scatter (e.g. 
Heavens & Jimenez 1999). One expects the TF scatter in bluer bands to be more sensitive 
to contributions of instantaneous star formation activity, while at (near-) infrared bands the 
scatter mainly reflects the convolved star formation histories. 

In order to verify that the basic trends in the scatter of the galaxy virial (TF and SL) 
relations are a manifestation of star formation histories, and thus persist at infrared bands, 
we need to expand our data base as we do below. 

2.3. The Infrared Tully-Fisher and Size-Luminosity Relations 

The advent of large-scale infrared surveys such as 2MASS provides us with reliable JHK 
luminosities, effective radii, and colors to construct (nearly) dust-free scaling relations. With 
a typical surface brightness limit of K ~ 20 mag arcsec" 2 , the 2MASS luminosity profiles are 
a full two magnitudes shallower than the typical J-band profiles in our samples or the SDSS 
brightness profiles. Yet, they yield scaling relations that are as tight as the ones derived at 
/-band. We now restrict our IR-extended analysis to the SCII sample for which the greatest 
range of luminosities and rotation speeds is found and thus for which the most reliable scaling 
relations can be derived. Because of the 2MASS magnitude limit of 13.5 K-m&g (e.g. Bell 
et al. 2003b), only the ~ 400 brightest SCII galaxies have measured Kron magnitudes. 

Figs. 11 to 13 show the 2MASS J-band scaling relations for SCII galaxies plotted against 
R e s, £ e ff, and J — K colors 4 . In the if -band one is notoriously more sensitive to airglow 
fluctuations (e.g., Jarrett et al. 2003), making this photometric band less attractive for TF 
studies. The if-band TFR is very similar to that in the J-band, but since the latter yielded 



4 Concentration indices from the 2MASS data base are known to be pathological (Bell et al. 2003b), as 
we verified ourselves. Consequently, we omit scaling relations with CI as a discriminant. 
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the smallest TF scatter, we focus on the J-band in what follows. The scatter in the SL and 
SV relations is significantly reduced compared to that in the /-band (cf. Table 2). Note, 
however, that the size parameter in the J-band is an effective radius rather than a disk 
scale length (which was not available from the 2MASS data release). Therefore, it is unclear 
whether this reduction in scatter reflects the usage of a more robust and rigorous galaxy size 
measure, or whether it is related to stellar population differences. 

As reported in Table 2, the (SCII) / and J TFRs have nearly the same log-slopes. For 
the SCII sample with 2MASS data, the infrared scaling relations are: 

V ohs oc L°/°, R d oc Lj 41 , .RdOcKbs 6 - (4) 

While the log-slopes for our SL and SV relations in Figs. 11 to 13 were computed with 
respect to R e s, they can be compared directly to i?<j since, for a pure exponential disk, 
R e g = 1.678.Rd- However, the J-band SL and SV log-slopes are significantly steeper than 
at /, possibly resulting from biases in the measurement of disk scale lengths versus effective 
radii. The former is sensitive to the disk fitting baseline while the other depends on the 
strength of the bulge-to-disk ratio. In order to verify this, we would need both types of radii 
measurements (either at J-band or in the infrared) which are not currently available; the 
assembly of such a data base is under way. Low and intermediate- redshift (z ~ 1) studies 
of SV relations, based on data from Courteau (1996; C97) and Vogt et al. (1997) and as 
reported in Mao, Mo, & White (1998), have SV log-slopes of 1.07 and 0.96 respectively, in 
good agreement with our /-band measurements. Dynamical considerations (§4) also suggest 
a log-slope of unity for the SV relation. More work needs to be done to ascertain the cause 
of the steep infrared size-dependent scaling relations. 

Figs. 11-13 confirm that the TF scatter is dominated by color and nearly independent 
of size and surface brightness. Previous reports of TF residual correlations with surface 
brightness may have generated confusion. Willick (1999) reported a correlation of /-band 
TF residuals with surface brightness and compactness for his "LP10K" survey of distant clus- 
ter galaxies using a moments fitting method to determine the exponential disk parameters 
even in the presence of irregularities in the galaxy light profiles. In an attempt to allevi- 
ate subjective fitting boundaries, Willick's computation used the entire surface brightness 
profile, including the bulge, and thereby biasing scale length and central surface brightness 
measurements. Applying his procedure, we can reproduce the putative surface brightness 
dependence of the TFR while proper bulge-to-disk (B/D) fitting techniques (MCH03) find 
none 5 . 



5 It is unfortunate that J. Willick is no longer with us to defend his approach. 
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The absence of surface brightness dependence of the TFR scatter has been verified by 
Sprayberry et al. (1995), Zwaan et al. (1995), and V01, and suggests a near-constancy of 
near-infrared (NIR) M^/L ratio for high surface brightness (HSB) and low surface brightness 
(LSB) galaxies of comparable total luminosities. In other words, if the surface brightness 
decreases, so do V b s and L in such a way that the TFR is independent of surface brightness. 
Thus, HSB and LSB galaxies belong to the same TFR, albeit with LSB galaxies showing a 
greater spread at lower luminosities. Any TF offset between low and high surface brightness 
noticed in the 5-band can be explained by mass-to-light ratio differences between the relevant 
stellar populations and, according to V01, the offset disappears at K. 

The SLR shows a very strong dependence on E off (as expected since S cff oc L/R^), and 
a weaker dependence on color. We see from Fig. 13 that bluer galaxies are on average larger 
than the mean, at a given luminosity, but we cannot state with any confidence that redder 
galaxies have smaller characteristic sizes, since these scatter evenly about the mean line. 

3. Residual Correlations of Scaling Relations 

The previous section identified color and surface brightness as the primary source of 
scatter in the TF and SL relations, respectively. In the spirit of CR99, we now examine 
correlations of the residuals from the mean virial relations in each sample. 

Following §2.1, we define the residuals for each object i as Ay{%) = y(i) —yf it (i). Fig. 14 
shows residual correlations for combinations of dlogV, dlogR, and <91ogL for each sample 
based on the /-band catalogs. The colored types have the same morphological dependence 
as in Figs. 1-4. The correlation residuals for the SCII sample with 2MASS J-band data are 
presented in Fig. 15; the point types are a function of J — K color. As in CR99, we apply 
a robust non-parametric test by rotating the set of residuals (Ayi(i), Ay 2 (i)) by various 
angles 9 to get (Ayi(i), Ay 2 {i)) and then applying a Spearman rank test (Press et al. 1992, 
§14.6) for correlations between the quantities (Ayi(i), Ay 2 {i)). The acceptable range of 
d(Ayi)/ d(Ay 2 ) can be calculated from the range of angles 9 for which the (Ayi(i), Ay 2 (i)) are 
not significantly correlated. The correlation slopes and their associated errors are reported 
in Table 3 and are shown at the bottom of each panel in Fig. 14. The linear correlation 
coefficients, r, which provides another means of assessing the strength of the correlation 
between each residual, are also shown at the top left corner of each panel. The slopes and 
correlation coefficients shown at the bottom of Fig. 15 are representative of the full data set 
(including all colors). 

There is close agreement in the correlated residual solutions of all samples, and including 
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the 2MASS SCII data set, with the general solution: 



d log V{L)/ d log R{L) 
d\ogR{V) / d\ogL{V) 
d log V{R)/ d log L{R) 



-0.07 ±0.05 



0.12 ±0.03 



0.29 ±0.02. 



These slopes are an average of the values for Shellflow, MAT, and the / and J-band SCII 
samples. Besides UMa, the slopes for the TF/SL residuals for all samples are negative 
(anti-correlated) and low. For the three large samples and at all colors, the TF/SL residual 
correlation is statistically different from d\ogV / dlogR = 0, though only weakly. Further- 
more, redder galaxies (top left corner of Fig. 15) lie above the null line whereas blue galaxies 
(bottom left) lie below. Note also that galaxies of all morphological types (and barredness, 
not shown here), scatter normally about the zero line. FAR00 have also computed the TF/SL 
residuals for the UMa sample (V01; see lower left panel of Fig. 16) and find no correlation 
for both for LSB and HSB galaxies, as do we for that sample (see Table 3). That result is 
however heavily weighted by the small size of the sample. 

The weak and strong surface brightness dependences of the TF and SL relations, respec- 
tively, have been used to infer the ratio of baryonic to dark mass in the inner parts of galaxies 
(CR99). If the gradient of NIR mass-to-light ratios is self-similar among galaxies of different 
brightnesses, CR99 conjectured that maximal disks 6 should be offset from the mean TF and 
SL relations in such a way that the fractional deviations, A \ogV(L) and A log-R(L), from the 
mean relations, V(L) and R(L), should be strongly correlated. The case of d log V / d log R = 
—0.5 would be expected for pure (maximal) exponential stellar disks (i.e., without a dark 
matter halo). Uncorrelated TF/SL residuals might then suggest that disk galaxies (of all col- 
ors, types, and barredness) are, on average, sub-maximal 7 and/or that a fine-tuning of halo 
and disk parameters washes out any surface brightness correlation in the TFR. Using the 
samples of bright late-type spiral galaxies of Courteau (1992) and MAT, with available R and 
/-band luminosity profile and Ha rotation measures, CR99 found a tentative correlation be- 
tween the TF residuals (measured in velocity space) and the SLR residuals. The TF/SL 
fractional deviations were found to be weakly correlated with d log V(L) / d log R(L) = 
-0.19 ± 0.05. With our expanded study, we now find d log V {L) / d log R(L) = -0.07 ± 
0.05, suggesting that the disk mass plays only a small role in setting the total velocity 



6 A stellar disk is "maximal" if it contributes more than 75% of the total rotational support of the galaxy 
at R = 2.2i?d, the radius of maximum disk circular speed (Sackett 1997). i?d is the scale length of the disk. 

7 The most massive spiral galaxies with V max > 200 kms" 1 may harbor disks with high mass fractions, 
in agreement with the maximum disk hypothesis (Athanassoula, Bosma, & Papaioannou 1987; Kranz, Slyz, 



& Rix 2003). 
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at 2.2R d . According to the galaxy contraction models of CR99 8 , the residual correlation 
d log V (L) / d log R(L) < -0.07 ± 0.05 corresponds to V disk /V total < 0.55 ± 0.10, indicating 
that spiral galaxies, independent of barredness or type, would be dominated by dark mat- 
ter with Md ar k/M tota i > 0.70 within 2.2R&. This result agrees with analytical models of 
cosmological structure formation (e.g. Dalcanton et al. 1997; MMW98; NS00; vdBOO; Cole 
et al. 2000), whereby dark halos dominate the potential even in the inner parts of galaxies, 
and with independent measurements of galaxy dynamics (e.g. Fuchs 2001; Trott & Webster 
2002; Kranz, Slyz, & Rix 2003). This result does not preclude the possibility of a stronger 
disk at radii smaller than 2.2R d . 

With dlogR(V) / dlogL(V) = 0.12 ±0.03, the SL residuals (middle column of Figs. 14- 
15) are also weakly correlated. A positive correlation may be real but not of fundamental 
interest. Consider a model where V(R) versus L(R) is a strict relation with no scatter, 
while R is a scattered quantity about the relation R(L). Then the scatter in both V(R) and 
L(R) comes solely from the scatter in R, and of course the residuals will be correlated, both 
determined by the value of R for a specific galaxy. Thus it is the TFR between V and L that 
gives rise to the V(R)-L(R) correlation. In a similar way, the weak SL correlation between 
R and L gives rise to a (weak) correlation between the residuals of R(V) and L(V). Models 
to explain the numerical value of the correlation residuals are presented in Paper II. 

The VL residual distributions (third column of Figs. 14-15) are tightly correlated with 
r ~ 0.75 (a little less at J) and relatively small slope errors; indeed, this relation is simply 
the TFR recast in differential form. 

Figs. 16 and 17 are expanded versions of the dlogV(L) / d log R(L) relations for the 
MAT, SCII, Shellflow, and UMa samples with J-band data, and the SCII sample with J-band 
data, respectively. Again, the correlated residuals are fully inconsistent with the prediction 
for a pure stellar exponential disk (short-dashed line), dlogV(L) / d log R(L) = —0.5 at 
Rd (CR99). The solid and long-dashed lines show the best fits from our Spearman rank 
test. The lower panels in Figs. 16 and 17 show TF and SL residuals with color. Trends 
of dlogV and dlogR with color residuals are shown in the right lower and central panels; 
the residual fits for the SCII+2MASS data are shown as solid and long-dashed lines in 
Fig. 17. Velocity residuals correlate tightly with color whereas size residuals are noisier and 
statistically consistent with zero. These trends match those seen in CR99 for the Courteau- 
Faber and MAT samples and ought to be reproduced by evolutionary models of disk galaxies 
(e.g. FAR00; vdBOO). If we account for the observed trends with color, and the simplest 
assumptions about single stellar population models, we find that the mean observed value of 



8 Our models have been independently confirmed by Bullock (2002; priv. comm.) 
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dlogV(L) I d log R(L) < —0.05 ± 0.05 yields an upper limit to the contribution of disk stars 
to the rotation curve at R = 2.2R d (CR99). 

4. On the Origin of the Disk Galaxy Scaling Relations 

The linear regressions of V^ bs on L and R d on L, mostly for the /-band results, yield 
the following average TF, SL, and SV relations: 

KbsOcL - 31 , /? d ocL a33 , tfdOcKbs - (5) 

Within these relations, at a given luminosity, we found that redder spirals tend to rotate 
faster and show hints of being more compact than bluer spirals. 

In the spirit of the spherical-collapse model (Gunn & Gott 1972) one can define the 
virial radius of a collapsed, virialized gravitational body as the radius inside of which the 
average density is a factor A times the average background density (the value of A depends 
on redshift and cosmology; see e.g., Eke, Cole & Frenk 1996; Bryan & Norman 1998). It 
is straightforward to show (see Appendix A) that with this definition the virial quantities 
R vir , M v i r (defined as the total mass within -R V i r ), and V v - U (defined as the circular velocity 
at i?vir) are related according to M vir oc V? T and i? vir oc V vir . Therefore, as long as M vir /L, 
Rvu/Rd, an d Kir/Kbs are roughly constant (i.e., are not significantly correlated with either 
of the virial quantities), these virial relations are all that is required to explain the observed 
scaling relations for disk galaxies. 

To realize the implications of this in terms of galaxy formation theory, we write 

M vir /L = (6) 

/bar £gf 

Here, = M^/L is defined as the disk mass-to-light ratio, /bar is the baryonic mass fraction 
of the Universe (/bar = ^b/^o ~ 0.1 for a ACDM cosmology), and e g f is the galaxy formation 
efficiency that describes what fraction of the baryonic mass inside the halo ultimately ends 
up in the disk. The virial relation M vir oc thus turns into the observed TFR L oc V^ 3 bs if, 
and only if, 

is not correlated with Vo bs and reveals an amount of scatter in agreement with that observed 
in the TF relation. 

In the standard model for disk formation, set forth by Fall & Efstathiou (1980), Dalcan- 
ton et al. (1997) and MMW98, the relation between R d and i? vir derives from the acquisition 
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and conservation of (specific) angular momentum by the baryons, when cooling to form a 
centrifugally supported disk. The total angular momentum of a virialized system is conve- 
niently expressed by the dimensionless spin parameter 

A = Jvir (8) 

V2M vir R V1I V vil 

(Bullock et al. 2001). For an exponential disk embedded in a dark matter halo, the total 
angular momentum is given by 

J d = 2tt f ^ E(R)V c {R)R 2 dR = M d R d V vir f v (9) 
Jo 

with S(-R) and V C (R) the surface density and circular velocity of the disk, and 

fv = e u — - du (10) 

JO V vir 

(MMW98) a dimensionless number. If disk formation occured under conservation of specific 
angular momentum, as generally assumed, one therefore obtains that 



-Rvir fv 

(MMW98). Thus, for the virial relation _R vir oc V vir to turn into the observed SV relation 
Rd oc V^bg, one requires that 

is not correlated with V a bs and has a scatter around its mean value in agreement with that 
in the observed SV relation. 

Similarly, in order to relate the virial relation between i? vir and M vir to the observed SL 
relation, the quantity 

Cs, = ±(¥) ll:i (13) 



fv VT d , 

cannot be correlated with galaxy luminosity, and has to reveal an amount of scatter that 
matches that in the SL relation observed. 

Thus, the challenge for galaxy formation theories is to understand how to meet the 
constraints on C T f, C S v, and C S l derived here. That this is not a trivial matter becomes 
clear if one takes into consideration that (i) both e g f and T depend extremely sensitively on 
the efficiencies of cooling, star formation and feedback, and are thus expected to be strongly 
correlated with M vir and A (see e.g., van den Bosch 2002a), and (ii) that both V ohs /V vir and 
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fv depend, in a convoluted way, on A, e g f, and on halo concentration, which itself depends 
on halo mass (e.g., Navarro, Frenk & White 1997). 

An additional challenge for galaxy formation models is to understand the dependence of 
the scaling relations on disk color. For instance, reproducing that color is the prime driver of 
scatter in the TF relation, requires that Ctf is correlated with color. In hierarchical structure 
formation scenarios, small mass haloes form earlier than their more massive counterparts. 
Since gas cooling is very efficient in low mass haloes at high redshift, these models generally 
predict that fainter galaxies (living in low mass haloes) are older, and therefore redder, than 
their more luminous counterparts. This is opposite to the observed trend, and has long 
been known to be an important problem for galaxy formation in CDM comogonies (see Bell 
et al. 2003a and MacArthur et al. 2003b for recent discussions). Because of this problem, 
simple disk formation models within this CDM framework naturally predict that the slope 
of the TF should be steeper in bluer bands (e.g., van den Bosch 2002a), again opposite to 
what has been observed. For the more massive galaxies to be older/redder, one could adjust 
the feedback prescription in order to vary the amount of gas accreted by spiral galaxies as 
a function of time. However, the energetic requirements for this to occur may exceed the 
energy reservoir in the stellar winds and supernovae and would require the action of an AGN 
(Bell et al. 2003a). The challenge is to realize this such that the constraints on Ctf, Csv, 
and Csl are not violated. 

In summary, we can conclude that the observed scaling relations of disk galaxies are 
predominantly a straightforward result of the basic dynamics of halo formation in the stan- 
dard cosmological picture and the simple, standard scenario of disk formation. They do, 
however, impose some interesting constraints on galaxy formation models, which we address 
in more detail in Paper II. 

5. Discussion 

By assembling some of the most extensive existing databases of galaxy structural pa- 
rameters, and adding new infrared luminosity and effective radii from 2MASS, we have been 
able to confirm or infer the following major observational results: 

• The scaling relations of spiral galaxies, obtained by linear regressions of Kbs on L and 
R d on L, are 

t t r0.31±0.02 p r0.34±0.02 p t/1.00±0.09 

Vobs OC Lj , K d <X Lj , K d OC V ohs 

This implies that the surface brightness £ oc L/R\ oc R d oc oc mV 3 . These scaling 
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relations seem to reflect the simplest possible model for disk galaxy formation based 
on angular momentum conservation of the baryons inside virialized dark matter halos. 

• The residuals of V(L) in the TFR are correlated with color, and uncorrelated with ef- 
fective surface brightness. The residuals of R(L) are correlated with surface brightness 
and marginally correlated with color. 

• Early-type (redder) disk galaxies rotate faster than later (bluer) types at a given lu- 
minosity. The effective radii of blue galaxies are marginally larger than the mean at 
a given luminosity, while the sizes of red galaxies do not show a significant systematic 
deviation from the mean. 

• The TF/SL residuals are weakly anti-correlated with <91og^(L) / d log R(L) = — 0.07± 
0.05. This confirms the results of CR99, suggesting that most spiral disks are sub- 
maximal. 

• The SV/LV residuals are weakly correlated with d log R(V) / Slog L(V) = 0.12 ±0.03. 

While the observed scaling relations of disk galaxies appear to be a straightforward 
result, the lack of correlation between the residuals from these relations of V Q h s and R& at a 
given L coupled with their color dependence is more puzzling. In Paper II, we argue that the 
small scatter in the TFR is consistent with being a reflection of the natural scatter in halo 
concentration c, as measured in cosmological simulations. This scatter is small because it is 
independent of the scatter in halo spin. Other explanations for the small size of the TFR 
scatter have invoked star formation and feedback processes (Eisenstein & Loeb 1996; Silk 
1997; Heavens & Jimenez 1999; FAR00), the small dispersion in halo and/or disk collapse 
times (Heavens & Jimenez 1999; vdBOO), and dynamical response of the halo to the disk 
assembly (NS00). Whether these, and other interpretations for the small scatter in the 
TFR, are all fully independent is still a matter of contention that we revisit, along with the 
interpretation of the major observational results above, in Paper II. 
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A. A General Derivation of the Tully-Fisher relation 

Below, we develop the theory for the TFR (see also MMW98; NSOO), and determine the 
fundamental nature of the TFR. This is because two simple predictions for the TFR have 
been used in the past: one with a log slope of 3 (—7.5 in magnitudes), and the other with a 
log-slope of 4 (—10 in magnitudes). Here we compare both predictions and show that they 
are related to one another. This discussion replicates §4, but giving more detail. 

Let us consider a dark-matter halo of mass M vir reaching virial equilibrium at a time 
corresponding to the cosmological expansion factor a = (1 + z) -1 . We assume that the halo 
properties at that time are relevant for the properties of the stellar component that has 
formed in this halo. 

For a body that is dominated by gravitational forces and is approximately in steady- 
state, the virial radius i? vir is defined in the spirit of the spherical-collapse model (Gunn & 
Gott 1972) by a given density contrast A relative to the mean universal density at the time of 
collapse (where mass shells are infalling for the first time), namely by M vir /i?^ ir oc A ~ 3 . At 
early times, when Q m ~ 1, the relevant density contrast is A ~ 178, while for the standard 
ACDM cosmology (with = 0.7 and Q m = 0.3 today) it rises to A ~ 340 today (e.g. Eke 
et al. 1998; Peacock 1999). We ignore the weak redshift dependence of the A factor in the 
following expressions 9 . 

If disks are embedded in dark halos we can define the total mass of the system (dark 
plus baryonic matter) as 

M vir = (Al) 

Here G is the gravitational constant, and R V1I is defined as the radius inside of which the 
average density of the system is ~ 200 times the critical density of the Universe, which at a 
given z is given by 

Pcrit = -J^G- (A2) 

One can thus write 

RviI = WW(zj' (A3) 



9 The maximum correction is obtained at low rcdshifts. For example, in the range z — — 2 the change 
is roughly A oc a 1 / 2 , which implies that a in the following expressions should be replaced by A _1 / 3 a oc a 5 / 6 . 
This is a weak effect, which becomes even weaker at higher rcdshifts where the bulk of galaxy formation 
occurs. 
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Here V v - W is simply the circular velocity at R VW1 and H(z) is the Hubble constant at the 
virialization redshift z of the halo. Note that a oc t 2 / 3 oc H(z)~ 2 ^ 3 , for an Einstein-de Sitter 
cosmology relevant at high redshifts, where t is loosely referred to as the time of formation 
of the halo. The three virial quantities at a thus define a one-parameter family: M vir oc 
a 3 /V v i oc a" 3 i? 3 ir . 

We now define the disk mass-to-light ratio 10 

T d , f , (A4) 

and write 

M d = / bar e gf M vir . (A5) 

Here /bar is the baryonic mass fraction of the Universe (/bar = ^b/^o ~ 0.1 for a ACDM 
model), and e g f is the galaxy formation efficiency that describes what fraction of the baryonic 
mass inside the halo ultimately ends up in the disk. Depending on the efficiencies of cooling, 
star formation and feedback, this can vary anywhere from zero to unity (see e.g., van den 
Bosch 2002a) 

Combining equations (Al) to (A5) yields 

L -ij0)r2 Vl - (A6) 

This is a TFR with log-slope 3 (see also Dalcanton et al. 1997; White 1997; MMW98; van 
den Bosch 1998; Syer, Mao, & Mo 1999). 

In contrast, several authors in the past have predicted that the TFR should have a 
log-slope of 4 (e.g., Sargent et al. 1977; Aaronson, Huchra, & Mould 1979; Salucci, Frenk, & 
Persic 1993; Sprayberry et al. 1995; Zwaan et al. 1995; Stil & Israel 1998). Their argument 
goes as follows: Assume that the disk surface brightness is described by an exponential 
function of the form 

E(r) = S exp(-r/ J R d ), (A7) 
with i?d the stellar scale length of the disk. The total luminosity of the disk is 

L = 2nZ R 2 d . (A8) 

Assuming that one measures the rotation velocity V^s at a radius r = si? d , and that the gas 
and stars in the disk move on circular orbits, 

VL ^ (A9) 



10 



We use the subscripts "d" and "h" to refer to the disk and dark matter halo, respectively. 
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with M(r) the total mass within radius r. If we define the total mass-to-light ratio 

T(r) = M ^ - Md(r) ± Mh(r) (A10) 
we can combine equations (A8), (A9) and (A10) to obtain 

L = 27^ ^ff^) Kt (AU) 

This is a TFR with a log-slope of 4, and as long as one measures the rotation velocity at 
a constant number of scale lengths, s, the scatter is determined purely by the variation in 
£ Y 2 (r). This was emphasized by Zwaan et al. (1995), who, from the finding that both HSB 
and LSB spirals follow the same TFR, concluded that LSB galaxies must have much larger 
values of T 2 (r) than their high surface brightness counterparts. Note that this "conspiracy" 
does not apply for the TFR of Eq. (A6), for which the zero-point is independent of the disk's 
surface brightness. 

The TFRs in Eqs. (A6) and (All) thus predict a different log-slope. The cause of this 
apparent paradox lies in the fact that the surface brightness, S oc V b s \~ 2 (MMW98), and 
Eq. (All) reduces to Eq. (A6). Thus, the two theoretical predictions of the TFR, with 
log-slopes 3 and 4, are directly related, but only the derivation for Eq. (A6) is strictly based 
on dynamics and it is the most fundamental of the two. 



B. SDSS Colors for Galaxy Classification 

Accurate optical colors based on digital imaging is available only for the Shellflow (V/) 
and UMa (BVI) samples. The MAT and SCII galaxy magnitudes were imaged only at 
/-band. In CR99, we computed a B — I color term for MAT galaxies using RC3 B-band 
fluxes but RC3 magnitudes are inherently uncertain with Am ~ 0.2 mag. With the release of 
SDSS/DR1 data (York et al. 2000; Abazajian et al. 2003), we can extract g and i magnitudes 
(as proxy for V and / magnitudes) for galaxies that overlap with Shellflow and SCII galaxies. 
We do this for Shellflow galaxies in order to map the SDSS system to the Cousins system, 
upon which Shellflow magnitudes are based. A final transformation can then be applied to 
SCII magnitudes. We compare in Figs. 18 & 19 instrumental Petrosian magnitudes from 
the SDSS first data release against instrumental magnitudes for 31 Shellflow and 39 SCII 
galaxies. The comparisons avoid any differences due to extinction and k- correct ions. Besides 
a zero-point offset, the g and % Petrosian magnitudes scale linearly with the V and I Cousins 
magnitudes in Shellflow and SCII. We compute transformations for the SDSS Petrosian 
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magnitudes into Cousins magnitudes for the SCII galaxies that overlap with the SDSS/DR1 
database. The instrumental "V" magnitudes for SCII were corrected the same way as the 
/-band magnitudes but using an extinction dependence Ay = 3.24, instead of Aj = 1.96, 
and the Poggianti (1997) formulation for the U-band k-correction. 
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Table 1. Redshift-Distance Galaxy Surveys 



Sample 


N 


Gal. type 


Phot, bands 


Mag/Diam. limits 


Redshift limits 


Rot. measure 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


MAT 


845 


field 


I 


£eso > 1-7 


< 7000 kins" 1 


Hi 


sen 


521 


cluster 


I 


12 < mi < 17 


[5000-19,000] 


Ha , Hi 


Shellifow 


340 


field 


v,i 


m B < 14.5 


[4500,7000] 


Ha 


UMa 


44 


cluster 


B,R,I,K 


m z < 14.5 


18.6 Mpc 


Hi 



Table 2. Virial Slopes and Correlation Coefficients 





MAT 




sen 


Shellilow 


UMa 




l0gV o bs(£/) 

log R d (Lj) 
\ogR d (V ohs ) 


0.31 log Li 
0.31 log L/ 
0.93 log Kbs 


(0.91) 
(0.72) 
(0.58) 


0.31 log Lj (0.90) 
0.34 log Li (0.74) 
1.00 log Kbs (0.61) 


0.31 logLj (0.86) 
0.35 log Li (0.58) 
0.64 log Kbs (0.42) 


0.32 \ogU 
0.181ogL 7 
0.52 log Kbs 


(0.94) 
(0.54) 
(0.54) 


hgV ohs (Lj) 

log R cS (Lj) 
log-Reff(Kbs) 






0.30 log Lj(0.87) 
0.41 log Lj(0.83) 
1.36 log Kbs(0.68) 









Table 3. Residual Correlations 



MAT 



sen 



Shellflow 



UMa 



d log V(L!)/ d log R(Lr) 

d\ogR(V ohs )/d\ogL(V ohs ) 

d\ogV(R d )/d\ogL(R d ) 



-0.04 ±0.02 (-0.26) 
0.12 ±0.03 (±0.30) 
0.30 ±0.02 (±0.78) 



-0.07 ±0.04 (-0.16) 
0.12 ±0.07 (±0.33) 
0.29 ±0.02 (±0.74) 



-0.08 ± 0.04 (-0.16) 0.03 ± 0.16 (±0.11) 
0.15 ± 0.06 (±0.35) 0.02 ± 0.08 (±0.03) 
0.30 ± 0.02 (±0.83) 0.29 ± 0.06 (±0.86) 



d log V(Lj) I d log R{Lj) 

d\ogR(V ohs )/d\ogL(V ohs ) 

dlogV(R d )/dlogL{R d ) 



-0.08 ±0.07 (-0.15) 
0.11 ±0.03 (±0.27) 
0.26 ±0.03 (±0.61) 
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log (L [L ]) log (V [km s" 1 ]) 

Fig. 1. — Scaling relations for SCII galaxies. Line widths are measured from Ha rotation 
curves and Hi line widths and disk scale lengths are measured using the "marking the disk" 
technique (see text). Luminosities are computed from fully-corrected /-band magnitudes. 
Unlike Dale et al. the magnitudes do not include any correction for morphological type 
dependence. The solid lines correspond to the linear regressions to the SCII galaxies, based 
on our data— model minimization technique. The inset gives the slope of the fit and the 
scatter of the ln-normal distribution. 
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Fig. 2. — Scaling relations for MAT galaxies. Line widths are measured from Ha rotation 
curves and disk scale lengths are measured from bulge-to-disk decompositions. Luminosities 
are computed from fully-corrected /-band magnitudes. The dashed line is the corresponding 
fit for the SCII sample. A significant zero-point difference is seen for the respective TFR 
fits. 



-34- 



Fig. 3. — Scaling relations for Shellflow galaxies. Line widths are measured from Ha rotation 
curves and disk scale lengths are measured from bulge-to-disk decompositions. Luminosities 
are computed from fully-corrected /-band magnitudes. The dashed line is the corresponding 
fit for the SCII sample (Fig. 1). 
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Fig. 4. — Scaling relations for UMa galaxies. Line widths are measured from Hi synthesis 
rotation curves and disk scale lengths are from bulge-to-disk decompositions. Luminosities 
are computed from fully-corrected J-band magnitudes The dashed line is the corresponding 
fit for the SCII sample (Fig. 1). 
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Fig. 5. — Combined TF and SL relations for the four data sets (SCII, Shellflow, MAT, 
and UMa) using absolute magnitudes (as conventionally plotted by the original authors). 
Zero-pointing differences clearly exist between each sample. 
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Fig. 6. — /-band magnitude residuals from the TFR for Shellflow galaxies. 



Fig. 7. — /-band magnitude residuals from the TFR for SCII galaxies. 
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Fig. 8. — Magnitude residuals, for the SCII and Shellflow samples, as a function of central 
surface brightness (/i = — 2.51ogS ), disk scale length (Rd), concentration index (CI), red 
color V — I, and optical morphological type (see text). The strongest and tighest correlation 
with TF residual is against galaxy color. 
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Fig. 9. — Same as Fig. 8 but for the MAT and UMa samples and using a B — I color 
term. CI's values were missing for the MAT galaxies and bulge-to-disk ratios, from Byun 
& Freeman (1995), were used as substitute. The slight trend of TF residual with /xq for 
the MAT sample is due to improper modeling of the galaxy bulge. Galaxy color shows the 
tightest correlation with TF residuals. 



-40- 



Fig. 10. — Correlation of galaxy observables (color, rotation speed, concentration index, and 
disk scale length) with morphological type. The point types represented different samples 
(green open squares: SCII (with g — % colors); red filled circles: Shellflow (with V — I 
colors); blue filled triangles: UMa (with B — I colors); the MAT galaxies were omitted for 
clarity.) Color is found to be most tightly correlated with morphological type. The degree 
of concentration, on the other hand is a very poor indicator of Hubble type. 
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Fig. 11. — Scaling relations for a subset of the SCII galaxies (N = 400) with J-band 
luminosities and effective radii, R e g in kpc. As in Fig. 1, the luminosities are not corrected 
for any morphological offset. The point types are a function of R e g. We see no striking 
dependence of TF scatter with R e g, other than the smaller systems liying slightly above the 
mean TFR. The solid lines are our best fits to the data computed as in §2.1. 
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Fig. 12. — Same as Fig. 11, but with point types as a function of effective surface brightness. 
The TFR is clearly independent of surface brightness; the scatter of the SLR is however 
strongly dominated by it. 
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Fig. 13. — Same as Fig. 11, but with point types as a function of J — K color. The scatter 
of the TFR is dominated by color. The SLR shows a weaker, though noticeable, dependence 
on color. 
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Fig. 14. — Correlation residuals from the mean relation at fixed virial quantity for the 
Shellflow, SCII, MAT, and UMa samples. Symbols are as in Fig. 3. Galaxies of all shapes 
and types (barred and un-barred, early and later types) obey the same surface brightness 
independent scaling relations. The residual distribution for d log V / d log R of all three 
samples is flat for barred and un-barred galaxies of all Hubble types and luminosities. The 
middle figures show no correlation between L and R, and the figures on the right side recast 
the TFR in its differential form. 
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Fig. 15. — Correlation residuals for the SCII sample with 2MASS J-band luminosities and 
effective radii, separated by J — K color to show the differences from red and bluer. The 
inscriptions in the bottom windows show the mean Pearson coefficients and residual slopes 
for all the colors combined. Redder/bluer galaxies lie slightly above/below the null line in 
the TF/SL relations (left panels). 
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Fig. 16. — TF/SL correlation residuals at fixed virial quantity for the Shellfiow, SCII, and 
UMa samples. Symbols are as in Fig. 1. The solid and long dashed line in the upper panel 
show the best fit and the statistically acceptable range of these slopes (at the 95% confidence 
level), respectively. The best fit is given by <91og^(L) / dlog R(L) = — 0.07±0.03. The short 
dashed line in the upper panel shows the maximal disk solution [d\ogV(L) / d log R(L) = 
—0.5]. The residual dependences on color (bottom) use the Shellfiow V — I data, the SLOAN 
g and i magnitudes transformed according to Figs. 18 & 19 for a subset of 39 SCII galaxies 
(see Appendix A), and the B — I colors for the MAT and UMa samples scaled to match the 
Shellfiow V — I range. 
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Fig. 17. — Correlation residuals for the SCII sample with J-band photometry. The line types 
are as in Fig. 16. The best fit is given by dlogV(L) / d log R(L) = — 0.08±0.07. The residual 
dependences on J — K color are shown in the bottom panels with the same line types as 
above for the right lower panels. Those fits are given by dlogV (L)/d(J — K) = 0.09 ± 0.04 
and d log R(L)/d( J — K) = — 0.19±0.15. The point types in the top figure refer to different 
colors (same notation as in Fig. 13) and the point types in the four bottom panels refer to 
different surface brightness levels (as in Fig. 12). 
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Fig. 18. — Comparison of raw SLOAN (Petrosian) vs Shellflow (Cousins) magnitudes. This 
confirms that the two sets are derivable from one another, independent of a color term at 
least over the small Shellflow magnitude range. 
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Fig. 19. — Comparison of raw SLOAN (Petrosian) vs SCII (Cousins) magnitudes. The linear 
transformation of magnitudes holds over nearly 3 magnitudes, independent of a color term. 



